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In the first 3.5 years of operations Fermi detected several sources with daily fluxes brighter than F{E > 
100 MeV) > 10~^ ph cm~^ the threshold set by the Fermi Collaboration for issuing an Astronomer Tele- 
gram (ATel). We focus the attention on these flaring sources, most of which are blazars, known to be extremely 
variable over the whole electromagnetic spectrum, from radio to gamma-ray energies and we investigate the 
properties of the selected sample comparing them to the general characteristics of the Fermi catalog sources. 



1. INTRODUCTION 

Since the beginning of its operations the Large 
Area Telescope (LAT) on-board the Fermi satellite 
has guaranteed a continuous monitoring of the MeV- 
GeV sky. Notably, the usual LAT operation mode, 
the scanning sky-survey mode, allows the coverage of 
the full sky every two orbits (about 3 hours) , therefore 
offering a great opportunity for rapid flare detection 
and follow-up observations. Complement to this all- 
sky monitoring capability of variability and transients 
detection is the invaluable effort offered by the Flare 
Advocate duty activity (also known as Gamma-ray 
Sky Watcher, FA-GSW), a scientific service belong- 
ing to the LAT Instrument Science Operations aiming 
to provide a first quick-look inspection and daily re- 
view of the gamma-ray sky observed by the LAT, car- 
ried out with continuity all the year round through 
weekly shifts. Sources that are found to exceed the 
predefined flux daily scale of F{E > 100 MeV) > 
10~^ ph cm~^ (which is a typical value used for 
defining a source in fiare at gamma rays) are promptly 
reported to the scientific community, e.g. by means of 
ATels. 

In the period from August 4, 2008 to February 4, 
2012, 91 of them fulfilled the ATel criterion as they 
underwent at least one flaring state overcoming the 
predefined flux threshold. We collect all these flaring 
events and we note that among them, the overwhelm- 
ing of flaring objects is composed by blazars. Actu- 
ally, although they represent a relatively rare sub-class 
of Active Galactic Nuclei (AGNs, 10% of the entire 
AGN population), these objects constitute the largest 
known population of gamma-ray sources as well as 
the great majority of sources detected by the Fermi 
LAT. We focus on the general characteristic of the 
sample by itself, paying particular attention to the 
general properties of the AGNs pointed out by the 
FA-GSW activity and provide a comparison with the 
AGNs seen in flaring state by the previous gamma-ray 
instrument, the Energetic Gamma- Ray Experiment 
Telescope (EGRE T), on-board the Compt on Gamma- 
Ray Observatory [Thompson et al. |l993 . 



2. SAMPLE COMPOSITION 



The objects that exceeded the fixed flux threshold 
that have been reported in Fermi Atels belong to dif- 
ferent source classes. Following the classification that 
has been used in the second Fermi catalog of AGNs 



[2LAC, Ackermann et al. 2011 the objects of this 



sample divide in: 68 Flat Spectrum Radio Quasars 
(FSRQs), 14 BL Lacs, 2 AGN, 4 AGU\ 2 Narrow 
Line Seyfert I (NLSyl) and one Radio Galaxy. 

Among them, 83 have a measured redshift and 7 
have not been detected in the second Fermi source 
catalog [2FGL, Nolan et al. 2012 . Considering the 



usual BL Lac classification based on the position of 
the SED peak, the 14 sources are divided in 6 ISP, 6 



LSP and 2 HSP las defined in Ackermann et al. 2011 



Worth to mention is that nine sources of this sample 
are detected at TeV energies: Ap Lib, S5 0716-1-74, 
BL Lacertae, 3C 66A, PKS 2155-304, 3C 279, 4C 
-t-21.35, PKS 1510-089 and NGG 1275. 



3. Fermi LAT DATA ANALYSIS 



The Fermi LAT is a pair conversion telescope with 
large effective area (~ 8000 cm^ on axis for E > 1 GeV) 
and large field of view ('-^ 2.4 sr at 1 GeV). It is op- 
timized for gamma rays in the energy range from 
20 MeV up to energies beyond 300 GeV. Full details of 
the instrument and descriptions of the on-board and 
g round data processing are provided in [Atwood et al.| 



[ 2009 , and information regarding on-orbit calibra- 



tion procedures is given by Ackermann et al. 2012 



Thanks to its high sensitivity and almost uniform sky 



^AGU refers to sources without a good optical spectrum 
or without an optical spectrum at all, whereas AGN refers to 
sources that are not conflrmed blazars nor blazar candidates. 
For more details we address the reader to lAckermann et al. I 
[20TTI. 
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coverage it is an ideal tool for niultiwavelength mon- 
itoring, gamma-ray flares detection and follow-up ob- 
servations. 

For our analysis we select all "source" class events 
located within 10° of each source of interest for the 
time interval reported in the correspondent ATel (one 
day interval). When a time interval larger than a day 
is indicated in the ATel, we look at the daily Auto- 
mated Science Processing (ASP) light curve and ana- 
lyze the data related to the first daily peak showed by 
the ASP light curve. To limit the contamination from 
Earth-limb gamma rays produced by the cosmic-rays 
interaction with the upper atmosphere, data are re- 
stricted to a maximal zenith angle of 100° and time 
periods when the spacecraft rocking angle exceeded 
52° are excluded. For each flaring episode we first 
use the pointlike tool to derive the best estimation 
on the position of the flaring source. Subsequently 
the analysis is performed with the standard analysis 
tool gtlike, part of the Fermi ScienceTools software 
package (version 09-27-01) available from the Fermi 
Science Support Center (FS SC^). An unbinned ma x- 
imum likelihood technique [Mattox et al. 1996 is 



applied to events in the energy range from 100 MeV 
to 300 GeV in combination with the post-launch in- 
strument response functions (IRFs) P7SOURCE_V6 
to derive the spectral parameters. 
The model for each Region of Interest (Rol) in- 
cludes the diffuse Galactic foreground emission and 
the isotropic diffuse emission: the former by means of 
the template gal_2yearp7v6_v0_trim.fits and the lat- 
ter, which accounts for both photons and residual 
charged particle background, by means of the tem- 
plate iso_p7v6source.txt^. Sources reported in the 
2FGL catalog and located within 15° of the target 
sources are incorporated in the model of the Rol as 
well. For these individual LAT sources and the sources 
of interest themselves we assumed a power-law model. 
In the fit the parameters of sources located within 10° 
radius centered on the source of interest are left free 
while parameters of sources located within the 10°-15° 
annulus are fixed and the isotropic normalization as 
well. We perform a first fit to estimate reliable start- 
ing values for the parameters in the Rol and remove 
from the model sources with too low significance. For 
this aim we use the Test Statistic (TS) which provides 
an estimate of the significance of the detection for each 
gamma-ray source in Rol. The TS value is defined as 
twice the difference between the log-likelihood func- 
tion maximized by adjusting all the parameters of the 
model, with and without the source. Sources with TS 
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Figure 1: Distribution of sources in Galactic latitude: 
solid line (gray area) represents the overall distribution 
corresponding to the 2LAC sample, dashed line this 
sample (normalized counts on vertical axis). 



lower than 5 are removed from the Rol model. A sub- 
sequent minimization procedure is applied to the data 
to derive the final values for the fluxes and spectral in- 
dices. When the target source is not signiflcantly de- 
tected, i.e. the TS value is lower than 25^, flux upper 
limits at 95% confidence level are calculated. 



4. GENERAL CHARACTERISTICS 

All the recalculated values and the complete out- 
comes resulting from this analysis will be presented in 
an upcoming publication, in the following we briefly 
report on the general characteristics of the sample. 
In particular we compare the properties of our sam- 
ple with those of the whole sources detected by Fermi 
considering the 2FGL and 2LAC catalogs. 

The sources of this sample appear to be equally 
distributed in the northern and southern hemisphere 
as evidenced by Fig. [T] which shows the histogram of 
the Galactic latitude distribution for the objects of 
this study (dashed line) and the ones reported in the 
2LAC clean catalog (solid line, gray area). The only 
discrepancy is found at low latitudes {\b\ < 10°), due 
to the fact that the 2LAC clean sample is build upon 
the a priori criterion to select only sources at high 
Galactic latitude, i.e. |&| > 10°. The redshift distribu- 
tion of this sample (dashed line) is presented in Fig. [2] 
and follows well the redshift distribution of the 2LAC 
clean sample (solid line, gray area). Fig. |3] shows the 
flux distribution of our sample as the ratio between 
the flux calculated for the time interval reported in 
the ATel and the average value reported in the 2LAC. 
We flnd that, on average, the sources of this sample 



^Web address: http://fernii.gsfc.nasa.gov/ssc/ 
^These templates are publicly available at the web address: 

http : / /f ermi . gsfc.nasa.gov/ssc/data/access/lat/ 

BackgroundModels . html 



^The square root of TS in the case of two degree of freedom 
is distributed as x^i therefore TS=25 roughly corresponds to 
4.6 sigma. 
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Figure 2: Redshift distribution: solid line (gray area) 
represents the overall distribution corresponding to the 
FSRQ 2LAC clean sample, dashed line the FSRQ of this 
sample (normalized counts on vertical axis). 
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Figure 4: Difference between the photon index reported 
in the 2LAC (which is computed over a period of 2 years) 
and the index measured in the time interval reported in 
the ATel. FSRQ are plotted with solid line, BL Lacs with 
dashed line; here only sources belonging to the 2LAC 
clean sample are considered. 




FIJX_Atel/FllK_2FQL 



Figure 3; Ratio between the flux (above 100 MeV) 
measured in the time interval reported in the Atel and 
the mean 2FGL flux; in the inset are shown the two main 
different classes: FSRQ solid line, BL Lacs dashed line. 



have fluxes a factor 11 higher than their 2LAC aver- 
ages. Despite this, there are some exceptional cases 
for which we measure a flaring flux > 40 times its 
average value (as in the case of 2FGL J0532. 0—4826, 
2FGL J1153.2+4935 and 2FGL J1848.6+3241). The 
flux distributions of the two main classes that make 
up our sample are plotted in the inset of Fig. [3} FS- 
RQs are represented by the solid line while BL Lacs 
by the dashed line. These distributions do not indi- 
cate strong differences between the classes, although 
we note that the higher flux increments are displayed 
only by the former. The distribution of the photon in- 
dex derived for the flaring state is presented in Fig. [4] 
evidencing the two classes of blazars: BL Lacs are 
plotted with a dashed line, and FSRQs with a solid 
line. The histogram represents the difference between 
the average 2LAC index and the index measured dur- 
ing the flare. While BL Lac indices remain almost 
constant, FSRQs show a hint of hardening their spec- 
tra during the flaring state. 



5. EGRET COMPARISON 

Among the 271 pointlike sources detected i n gamma 
rays by EGRET [see e.g. Mukherjee 2001 66 were 



blazars plus a radio galaxy (Centaurus A), and some of 
them showed variability, being observed several times 
at different amplitudes Hartman et al. 1 1999 . From 



a comparison of the AGNs detected by EGRET with 
the objects observed in flaring activity by Fermi LAT 
during the first 3.5 years of operation, results that 29 
of the 67 EGRET sources have been detected in high 
gamma-ray state by Fermi. Focusing on the gamma- 
ray sources detected by EGRET with a flux F{E > 
100 MeV) > 10"*^ ph cm-2 s'^, 10 of the 16 flaring 



sources observed by EGRET [see e.g. Mukherjee [2001 



have displayed intense gamma-ray activity also during 
the Fermi operation. This is an indication that the 
high activity of these sources continued at least over 
two decades, even if in some cases the source under- 
went a quiescent phase before increasing again its ac- 
tivity. A typical example is PKS 0528-1-134, that was 
observed by EGRET during di fferent flaring ep isodes 
between 1991 and 1997 [Mukherjee et al. |1999| , while 
it remained quiescent during the Fermi observation 



up to a new flare in June 2011 D'Ammando 2011 
On the contrary, e.g. 3C 279 was active for most of 
the time during both the EGRET and Fermi observa- 
tions. However, most of the flaring AGNs detected by 
the LAT did not show significant gamma-ray activity 
during the EGRET observations. In this context, it is 
indicative the behavior observed on 3C 454.3 and PKS 
1510—089, the two brightest extragalactic blazars in 
the Fermi sky. During the EGRET observations only 
a moderate activity was observed from these sources; 
on the contrary they were characterized by an intense 
and dynamic behavior during the Fermi monitoring, 
reaching a remarkable gamma-ray fiux higher than 



F{E > 100 MeV) > 10"^ ph cm~^ s ~^ Sanchez and 
Escande ||2012[ [Hungwe et al. [2011 
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Indeed, taking into consideration the different oper- 
ating mode of EGRET and LAT, i.e. inertial pointing 
versus all-sky scanning survey mode, it cannot be ex- 
cluded that some strong gamma-ray flares have been 
missed for lack of observation coverage in the EGRET 
era. Notwithstanding, the finding that some blazars 
detected by EGRET and not by Fermi suggests a pos- 
sible significant change of activity in some objects, at 
least between the epochs of EGRET and Fermi oper- 
ation. 



6. CONCLUSIONS 

During the first 3.5 years of Fermi operations 91 
sources displayed intense gamma-ray activity with 
fluxes exceeding F{E > 100 MeV) > 10"^ cm'^ s'^. 
We studied this sample of sources comparing their 
characteristics to the 2FGL and 2LAC samples and 
we showed that their properties are similar to the 
main Fermi LAT samples. It is worth to notice that 
this sample constitutes only 10% of all the sources de- 
tected in 2 years of Fermi operations. Therefore, up 
to now, only a small percentage of the 2FGL sources 
displayed variability reaching extreme flux values, and 
those are basically all blazars. All the flaring episodes 
have been analyzed in details and the correspondent 
results will be presented in the upcoming publication 
where also further considerations will be addressed. 



Acknowledgments 

The Fermi LAT Collaboration acknowledges sup- 
port from a number of agencies and institutes for 
both development and the operation of the LAT 



as well as scientific data analysis. These include 
NASA and DOE in the United States, CEA/Irfu 
and IN2P3/CNRS in Prance, ASI and INFN in 
Italy, MEXT, KEK, and JAXA in Japan, and the 
K. A. Wallenberg Foundation, the Swedish Research 
Council and the National Space Board in Sweden. 
Additional support from INAF in Italy and CNES 
in France for science analysis during the operations 
phase is also gratefully acknowledged. 



References 

Atwood, W. B. et al., 2009 "The Astrophysical Jour- 
nal" 697,1071-1102 

Ackcrmann, M. ct al., 2012 "The Astrophysical Jour- 
nal Supplement" 203, 4 

Ackermann, M. et al., 2011 "The Astrophysical Jour- 
nal", 743, 171 

Nolan, P. L. et al., 2012 "The Astrophysical Journal 
Supplement", 199, 31 

Hartman, R. C. et al., 1999 "The Astrophysical Jour- 
nal Supplement", 123, 79-202 

Mukhcrjee, R., 2001 "American Institute of Physics 
Conference Series", 558, 324-337 

Mattox, J. R 1996 "The Astrophysical Journal" 461, 
396 

Mukherjee, R. at al., 1999 "The Astrophysical Jour- 
nal", 527, 132-142 

D'Ammando, F., 2011 "The Astronomer's Telegram", 
3412, 1 

Sanchez, D. and Escande L., 2012 "The Astronomer's 
Telegram", 3041,1 

Hungwe, A. et al., 2011 "The Astronomer's Tele- 
gram", 3694, 1 

Thompson, D. J., 1993 "The Astrophysical Journal 
Supplement" 86, 629-656 



eConf C121028 



